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ABSTTUCT 

The  function  E  =  C  +  iS  of  the  title  arises  in  diffraction  problems 
involving  gratings  and  in  allied  antenna,  problems.   We  present  curves  of 
S  vs.  C  as  fionctions  of  6  for  IT  =  1  to  10  and  N  =  10^,  lO"^,  10  .   It  is 
found  that  JE]    always  occurs  for  6  =  0,  and  the  curves  from  6  =  0  to  2tt 
show  N-1  loops;  for  large  N,  JEJ    occurs  for  6  =  tt.  Curves  of  C  and  S 
83  functions  of  6  are  also  included. 
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The  function  ' 

N 

(1)   l](N,5)  =C  n"^/^  e^""^   =  C  +  iS, 
n'=l 

where  C  and  S  are  real,  is  the  limiting  form  of  a  Schlbmilch  aeries  which 

arises  in  the  multiple  scattering  solution  of  the  grating  of  parallel 

cylinders  i"   It  also  arises  in  allied  antenna  problems  involving  linear 

arrays  of  coupled  elements.   In  these  problems,  S  is  the  phase  difference 

of  neighboring  elements  and  N  is  the  number  of  neighbors  taken  into  account. 

The  function  E,  which  involves  two  variables,  is  essentially  the  sum  analog 

of  the  Fresnel  integral. 

For  small  values  of  N  we  have  evaluated  S  by  graphical  methods.  Figures 
1  to  10  are  curves  of  S  vs.  C  for  N  "  1  to  10;  the  variable  5'  ranges  from  0 
to  180**.  These  curves  if  taken  from  6=0°  to  360°  have  N-1  loops.  With 
Increasing  N  the  train  of  loops  forms  a  "pigtail"  spiral,  and  JE|    tends  to 
occur  for  6  =  180°,   i-^lnjax  ^^*'®^^  occurs  for  6=0, 

Figures  11  to  20  are  crxrves  of  C  and  S  as  functions  of  6  for  N  =  1  to  10, 
The  curves  for  N  =  1  are  simply  the  cosine  and  sine  curves;  the  corresponding 
curves  for  N  >  1  have  essentially  the  same  symmetry  properties  as  the  cosine 
and  sine. 

Figures  21  to  24  are  vector  diagrams  for  several  constant  values  of  6  and 
N  increasing  from  1  to  10.  For  a  fixed  large  value  of  N,  graphs  of  this  type 
may  have  essentially  one  of  three  forms  as  determined  by  the  constant  6, 
For  6=0  (or  tt)  the  plot  is  a  straight  line  or,  more  generally,  for 
6   1(0+  tt/2)/N,  0  =  tan~^(S/C),  the  plot  is  a  'simple  arc'  and  the 
resultant  of  the  vectors  is  much  larger  than  xinity,  as  in  Figure  21,*  For 
C  t     +  Tr/2)/N  <  ^  i  tt/z,  the  plot  is  a  'spiral'  as  in  Figures  22  and  23. 


Cl3  7.  Twersky,  J,  Appl,  Phys,  23,  1099  (1952).  See  also  paper  by  S.Karp;  in 

press. 
*  Tne  transition  value  of  6  is  deduced  as  follows:   for  the  simple  arc  case 

to  obtain  we  require  that  the  N-th  vector  make  an  angle  less  than  90  with  the 

resultant,  i.e.,  if  the  angle  is  greater  than  90   the  graph  has  already  begun 

to  spiral;  hence,  for  large  N,  the  transition  from  simple  arc  to  spiral  occurs 

for  N  6  =  0  +  T7/2.   We  are  grateful  to  W,  Magnus  for  suggesting  this. 
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For  tt/2   <   5   <  tt,    the  plot  forms    a  Jag^d    'star  patch'    which  intersects 
itself  maay  times   aad  has   a  resultant   less    than  unity;    cf .   Fig.   24. 

For  large  values   of   N,   a  more   rapidly  convergent  representation  of   ^ 

suitaole  for  nxmerlcal  calculations  may  be  obtained  by  means   of   the    'Euler 

SI 
sum  formula"^      thus 

(2)        C(N.6)   =^/^   [C(N6)    -^(6)]-^(cos    6  ^  ^M  ) 

+     I_   (Cfis  J  +   6   sin   6   _  cos  JfS_  _        6  a^n  ^6   ^  ...    , 

+    i-   (BiTLi  _    5  ^^^  5  _  sin_N6_     ^     ^gos   N6  .   ^    ...    , 
^^  ^  2^^  ^^ 

where  (^  and  ^  are   the   tabulated  Fresnel  integrals  J -^J      These  are  rapidly 
convergent  for  N  large  or  6    small.     Thus   for  S    =  0,   we  have 

C   =  Zn~^'^^2>^  -  1.46  +  l/2<^  ;    for  N  >   5O,    this   expression  yields 

results   agreeing  to  at   least   three   places   with  those   obtained  by  adding 

tabulated  values   of     n"''"'^.     We  find  C  =  18,6,   6I.7,   198.5     for  N  =  10^, 

3         4 
10    ,    10    ,    respectively. 


For  valvies   of  6    near  tt  we   write     6    =  tt  -  c    and  rewrite    (l)   as 
(3)     £( 


E/    ^nQ  -ine  v-"       -lQ«         r^       -inc/2 

n=l  «/n  n=1       *  n  n— 1      \' n 


n=l  vn  n=l    *■  n  n-1    vn 

=  -E(N,-<)   +^'F"fi(N/2   -  €/2), 

or  as 

(3')      S(N,6)    =  -  C(N,€)   +  ^^  C(N/2.    ^/z)   +   i    [S(N,e)    -^2   S(N/2.    c/2)] 
in  order   to  use   the   representations   of    (2).      Thus  for   5  =  tt, 

|_2l      See,  for  example,   L,  Bieberbach,   Lehrbuch  der  Funictlonen   Theorie    (Chelsea 

Publisning  Co.,  N.Y. ,   1945)   Vol.   I,  p,   30?, 
[3i     ii.  Jahjoke  and  F.   Emde     ,  Tables   of  Functions    (Dover  Publications,   N,Y, , 

1945)  p.   35. 
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n=l 

2    3    ^ 
Firnres  25,  27,  29  are  plots  of  S  vs.  C  for  N  =  10  ,  10"^,  10  respectively 

snd  Tlrares  26,  23,  "Kl   are  the  corresponding  plots  of  C  and  S  as  functions  of 

u  =  (2N6/tt)-'-' ^.   With  increasing  N,  the  'axis'  of  the  spiral  tends  to  "become 

parallel  to  the  line  C  =  S.  Figures  31,  32  and  33  are  enlargements  of  the 

2    3^ 
toils  of  the  sT^ir-'ls  for  N  =  10  ,  10'  and  10  respectively.  Figure  JM   is 

the  analogous  CTorve  for  N  =  oo  o"btalned  "by  Sollfrey  and  Shmoys.-  -^      It  can 

"be  seen  that  in  this  range  of  6  (spv  6  =  kO°   to  180°)  the  function  is 

practically  Independent  of  N,  provided  that  N  is  large. 

The  values  of  the  graphs  are  accurate  to  about  five  percent. 

To  illustrate  the  use  of  the  fimctlon  E(N,6)  we  consider  the  scattering 
of  a  plane  electromagnetic  wave  normally  incident  on  an  infinite  grating 
of  parallel  cylinders  with  radius  a  and  spacing  "b.   If  we  assume  that  the 
incident  wave  is  polarized  parallel  to  the  cylinders,  and  that  Ka«l,  Kb»l 
(K  =  2tt/\),  we  may  write  the  total  field  (the  amplitude  of  the  electric 
intensity)  as 

W=-00 

For  perfectly  conducting  cylinders  we   impose  the  "boundary  conaition 
>i/  =   0  at   r     =  a  and  ohtain 

V 

00 


0  =  1  +  a[h^^^    (Ka)  +   2II    HjJ^^    (nKb)]; 


n=l 

we  have  employed  erpClKaccs  0)«   1,     Q  =  tan~    (y/x),  for  simplicity. 
Hence  we  may  write 

(5)     A  =  A^/[l-   2Aq  C  H^^^(nXb)]  .     A^  =  -l/H^^^(Ka),   Tto«l, 

where  A   is  essentially  the  scntterlng  coefficient  of  an  isolated  cylinder. 


The  scattered  component  of  (^)  is  then 


[5] 


[U]      W,  Sollfrey  and  J.  Shmoys,  Diffraction  of  Electromagnetic  Waves  "by  a 
PlFne  Wire  Grating.  EM-^2  (Math.  Res.(Jrp.,  N.Y.U. ,  N.Y.,  1952). 

r?}   Solutions  of  this  form  were  obtained  previously  by  R.  Honerjager, 
Ann.  Physik  k,   25  (19^8)  snd  W.  Frnnz,  Z.  ang.  Phys .  £,  ^16  (19^9). 
The  present  derivation  is  due  to  S.  Karp  and  J.B.  Keller. 


-  k  - 


The  numerator  m;  ,   is  the  result  obtained  if  each  cylinder  is  treated  as 

completely  isolated  from  the  others,  i.e.,  the  'single  scattering'  solution. 
The  denominator  D  may  be  expanded  as  a  power  series  and  ^  readily  interpreted 
physically  in  terms  of  hii^her  order  scattering  processes;  yy   is  thus  an 
approximate  'multiple  scattering'  solution  for  the  problem, 

^ince  Kb>>l  we  may  employ  the  asymptotic  form  of  H.   (nKb)  to  obtain 

(7)   M/  /u/'  =  D"^  =  [1-  2A^(2/TTKb)^/2  e'^"^^  E(N,6)]"^ 

where  N  is  the  number  of  neighbors  we  take  into  account,  and  6  =  Kb  -  Zm^r 
is  the  phase  difference  between  waves  from  successive  neighbors.  Graphs 
of  i^l~^     and  additional  illustrations  will  be  found  in  reference  £1]. 


71p-ure8   1   to   10 

GraT)h3  of  S  vs.   C  for  TT  =  1   to  10;   the  parn::ieter  6  ranges  from  0     to 
180°,    the  ciirvea  "beln^  symmetrical  to   the  real  axis. 
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Figure  2 
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J'igure  5 


I^igTire  6 
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Figure  9 


Figure   10 
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TifTMTeQ   11   to   20 
C  and  S  as   functions   of   6  for  N  =   1  to   10 
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Figure  12 
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Figure  13 


Figruxe   1^ 
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Figure  18 
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Figure   20 
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Tir-areQ    21  to  2^ 

Vector  plots  for  6  constant  and  N  increasing 
from  1  to  10 


Figure  21 
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figure  22 
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Figure  25.      C  vs.   S  fts  fxmctions  of  u  for  IT  =  10 
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l^igTire   26.      C  find  S  as   functions   of  u  for  K  =  10 
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Flgrure   27.      C  vs.   S  as  functions   of  u  for  N  =  10"^ 
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"FlCTire   28.      C  and   S  as   functions   of  u  for  N  =  10"^ 
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Fl^re   29.      C  vs.    S  as  fimctlons   of  u  for  H  =  10 
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yig^xre   30.      C  anc^    S  as   fainctlons   of  u  for  IT  =  10 
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Y^fn^TP.a    "^1   to   3^ 
Enlargements   of  the   C  vs.    S  c^ja-res   for  the  vicinity 

of     6    =    TT 
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Figure  3^ 
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Graphs  of  the  function... 


N.Y.U.  Courant  Institute  of 
Mathematical  Sciences 

4  Washington  Place 
New  York  3,  N.  Y. 


